Entanglement spectra of critical and near-critical systems in one dimension 
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The entanglement spectrum of a pure state of a bipartite system is the full set of eigenvalues of the reduced 
density matrix obtained from tracing out one part. Such spectra are known in several cases to contain important 
information beyond that in the entanglement entropy. This paper studies the entanglement spectrum for a variety 
of critical and near-critical quantum lattice models in one dimension, chiefly by the iTEBD numerical method, 
which enables both integrable and non-integrable models to be studied. We find that the distribution of eigen- 
values in the entanglement spectra agrees with an approximate result derived by Calabrese and Lefevre Q|] to an 
accuracy of a few percent for all models studied. This result applies whether the correlation length is intrinsic or 
enerated by the finite matrix size accessible in iTEBD. For the transverse Ising model, the known exact results 
for the entanglement spectrum are used to confirm the validity of the iTEBD approach. For more general 
models, no exact result is available but the iTEBD results directly test the hypothesis that all moments of the 
reduced density matrix are determined by a single parameter. 

PACS numbers: 64.70.Tg, 03.67.Mn,75.10.Pq 



I. INTRODUCTION 

The past decade has seen a fruitful interplay between prob- 
lems in condensed matter physics and concepts from quantum 
information theory. Entanglement entropy, as reviewed below, 
has been studied for a wide variety of quantum many-particle 
systems and has given some considerable insight into both 
quantum criticality and quantum topological phases. More 
recently it has become clear that for some purposes the en- 
tanglement entropy is an insufficient characterization of the 
quantum information in a quantum state. Both for topological 
phases and for quantum critical systems, several recent stud- 
ies have shown that more information about entanglement, the 
full "entanglement spectrum", is required to describe impor- 
tant physical properties. The purpose of this article is to exam- 
ine entanglement spectra at and near several quantum critical 
points in one dimension, chiefly by numerical methods. This 
serves as a test of existing theoretical predictions. The main 
subject of this paper is a set of numerical tests of the validity 
of the approximate distribution of eigenvalues in the entangle- 
ment spectra derived by Calabrese and Lefevre JH]. 

The entanglement entropy of a pure state |W) of a bipartite 
system AB is defined as the von Neumann entropy of reduced 
density matrix pa, obtained by tracing out subsystem B. The 
result is symmetric with respect to A and B: 



-Tr pa log pa = -Tr Pb logps = S B . 



(1) 



This symmetry follows from the Schmidt decomposition, 
which is convenient for our later discussion. The Schmidt 
decomposition yfl for a bipartite system expresses the orig- 
inal wavefunction as a sum of product states of wavefunc- 
tions for the two parts of the system, with the bases chosen 
for those two parts so that the terms in the sum involve or- 
thogonal wavefunctions: 



|*) = ^A n |$ IlA )|$„B> 



(2) 



spaces of the subsystems A and B, and the A„ are positive. 
The Schmidt decomposition contains more than one term for 
entangled states, and the entanglement entropy is 



S 



uj n logw„ 



(3) 



The l$ n A) and form orthonormal bases for the Hilbert 



Note that the eigenvalues of the reduced density matrix are the 
squares of the coefficients appearing in the Schmidt decompo- 
sition and uj n = A 2 . 

Clearly the entanglement entropy is a particular combina- 
tion of the eigenvalues of the reduced density matrix. The "en- 
tanglement spectrum", i.e., the full set of eigenvalues of the 
reduced density matrix, clearly contains additional informa- 
tion in principle. For topological phases arising in the quan- 
tum Hall effect, the entanglement spectrum has been argued to 
contain information about topological order 0,|5t] beyond that 
of the entropy, which contains information about the "quan- 
tum dimensions" of the topological phase 10,0]. The entropy 
has been a powerful tool in numerical studies to identify topo- 
logical phases Jl,|^]. 

The value of entanglement spectrum for critical points, as 
studied in this paper, is slightly different from its use in mea- 
suring topological order. It has been argued by Calabrese and 
Lefevre QJLD that quantum critical points have an entanglement 
spectrum with a universal one-parameter dependence; since 
this one parameter can be taken to be the entanglement en- 
tropy, one might think that for quantum critical points, the 
entanglement spectrum contains no more useful information 
than the entanglement entropy. However, the detailed form 
of the entanglement spectrum is important; as an example, 
the Calabrese-Lefevre approximate form of the entanglement 
spectrum was used recently as the basis of a theory for "finite- 
entanglement scaling", i.e., how the finite entanglement avail- 
able in a matrix product state representation affects the ob- 
served behavior at quantum critical points. For this applica- 
tion, understanding not just the entanglement entropy but the 
full spectrum is essential, and the nature of this spectrum near 
a one-dimensional quantum critical point is the subject of this 
paper. 
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The focus on one-dimensional quantum critical points re- 
sults from two technical considerations, one numerical and 
one analytical. Numerically, the matrix-product-state class of 
numerical methods (DMRG and its descendants), reviewed 
in the following section, are especially powerful in one di- 
mension and capable of giving essentially exact ground-state 
properties for many systems away from criticality. Analyti- 
cally, conformal invariance of the d + 1-dimensional path in- 
tegral applies at most interesting critical points in any spatial 
dimension d, but only for d = 1 is the (local) conformal al- 
gebra infinite-dimensional. This infinite-dimensional algebra 
closely constrains the behavior of one-dimensional systems 
close to critical points and has led to numerous predictions 
about entanglement entropy at such critical points, many of 
which have been confirmed numerically. 

In the following section, we review the relevant analytical 
results on entanglement entropy and entanglement spectrum 
at one-dimensional critical points. Section III briefly reviews 
matrix product states and the Infinite Time-Evolved Block 
Decimation fioll (iTEBD) and presents our numerical results 
on several quantum critical points, including some integrable 
cases as a check and also some non-integrable cases. 

II. ENTANGLEMENT ENTROPY AND ENTANGLEMENT 
SPECTRUM 

An important early example of overlap between condensed 
matter physics and quantum information was the realization 
that "most" quantum critical points in one dimension (those 
with conformal invariance) have universal entanglement en- 
tropy determined by the central charge of the conformal field 
theory (CFT). Since central charge was already known to 
determine properties of the critical point related to entropy, 
such as the free energy at small nonzero temperature, it is in- 
tuitively plausible that the entanglement entropy is also re- 
lated to this quantity. For a partition of an infinite one- 
dimensional chain into a contiguous region of N sites and the 
infinite remainder, the entanglement entropy scales for large 

n as mi miii 

SW~|logiV. (4) 

When the system is away from criticality, this logarithmic di- 
vergence is cut off by the correlation length £, measured in 
units of the lattice spacing a: 

SW->|]og(f/a). (5) 

These entanglement entropies are universal since they depend 
only on central charge; different critical points in the same 
universality class have the same entanglement entropy. An el- 
egant derivation of these results lfl3ll is via a replica trick that 
obtains the entanglement entropy using the theory defined on 
a particular Riemann surface; the central charge then appears 
via the mapping of an ordinary plane to this Riemann surface. 
It is believed that the full entanglement spectrum that leads 
to this entropy is not universal but nearly so y]], as now ex- 
plained. Since the entanglement spectrum is not universal, it 



should not follow from the conformal mapping method that 
makes clear the universality of the entropy. 

There are only a few results on the behavior of entangle- 
ment entropy at translation-invariant criti cal p oints above one 
spatial dimension H [H H [H, M, M, El EB and in 
general the entanglement spectrum has not been considered. 
Hence we focus on the entanglement spectrum at or near a 
one-dimensional quantum critical point with conformal in- 
variance. The entanglement spectrum can be obtained from 
traces of all powers of the reduced density matrix: 

R a = Trp% (6) 

In order to make the entanglement spectrum well-defined, we 
have to introduce a length scale L e s coming from either a 
block size or a correlation length, just as for the entanglement 
entropy. When this length scale is much larger than any mi- 
croscopic scale, a scaling analysis jTl] predicts 

n _ r -c(a-l/a)/6 , 7 s 

where c a are nonuniversal constants. Following Ref. []1 we 
introduce the parameter b and rewrite this as 

R a = c a e- b ^- 1 ' a \ 6=^logi off (8) 

6 

and can obtain the form of the entanglement spectrum imme- 
diately under the assumption that the c a are constant. This 
assumption of constant c a was shown to be a good approx- 
imation for the XX model and even on the whole XX Z 
line (more details are given later). However, it seems valu- 
able to evaluate the validity of this assumption more gener- 
ally, including for models that are not integrable. Note that 
the parameter b is simply related to the entanglement entropy 
S obtained for this distribution: b = 5/2. For conformally 
invariant models, b is related to the largest eigenvalue of the 
density matrix b = — logw max . This is the so called single- 
copy entanglement EH EH E3l • We found numerically, that 
this relation holds in good approximation if we slightly detune 
the system from criticality. 

Under the assumption of constant c a , a Laplace transform 
leads to a universal one-parameter distribution of density ma- 
trix eigenvalues near the critical point Qj]. The parameter ap- 
pearing in the distribution, which we take to be b, combines 
both an "intrinsic" property of the critical point (the central 
charge) and an "extrinsic" one (either the correlation length 
or the block size, whichever is smaller). An intuitive way to 
express this distribution is via the mean number of eigenvalues 
larger than a given value u, 

n(w)= / dwP{w) = /o(2V61og(w max /a;)). (9) 

J u 

Here the largest eigenvalue w max is determined by the one 
parameter b: b = — log w max . As b increases, the tail of the 
distribution becomes longer. Examples of this distribution are 
shown in FIGS.|2]-|4]below. The entanglement spectrum of ID 
systems away from criticality has also been discussed in the 
context of majorization theory in Ref. EH 
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III. NUMERICAL METHODS AND RESULTS 
A. Methods 

Our main results are obtained using the Infinite Time- 
Evolving Block Decimation (iTEBD) algorithm [1 1 Ofl to study 
several one-dimensional Hamiltonians with translationally in- 
variant ground states. This algorithm can be viewed as 
a descendant of the density matrix renormalization group 
(DMRG) algorithm JUt Both DMRG and iTEBD construct 
trial wavefunctions that are "matrix product states" (MPS) 
l27tl . We review the basics of such states in order to under- 
stand their connection to the entanglement spectrum. 

We consider a system with N sites and periodic boundary 
conditions, where each site has d orthogonal states. Any pure 
state of the system is a superposition of the product basis states 
\s1S2 ■ ■ ■ sn) = \{s}) where 1 < s; < d. A MPS for such a 
system has the form 



E 



Tr 



Sl,...,Siv=l 



4 1 ' 



A [N] 
s N 



\si)...\a N ). (10) 



For each site i there are d matrices A l s . of a finite dimension 
X x x- A product (unentangled) wavefunction for a chain of 
particles/spins is obtained by multiplying together scalar am- 
plitudes for the particle/spin state at each site. The MPS gen- 
erates entanglement by using matrices for each site and can 
store an increasing amount of information as the matrix di- 
mension x increases. For fixed x, there are at most x nonzero 
eigenvalues of the reduced density matrix, and hence the max- 
imum possible entanglement is log x- 

The advantage of requiring translational invariance is that 
there are only finitely many different matrices A\.. The 
iTEBD algorithm finds an approximation of the ground state 
by performing an imaginary time evolution of a randomly 
chosen initial MPS iV'init) w ith fixed dimension x- It is 
therefore required that |V'in; t ) h as a non-zero overlap with the 
ground state. Because this method always constructs wave- 
functions for the infinite system, its errors result from finite 
entanglement rather than finite size, and one application of 
the entanglement spectrum was to develop a theory 12811 for 
the "finite-entanglement scaling" of resulting errors 12911 : the 
effective correlation length induced by finite entanglement for 
a critical Hamiltonian scales as £ oc x K , where the exponent k 
is determined by the central charge at the critical point. For the 
iTEBD simulations presented in this paper, we used a second 
order Trotter decomposition for the imaginary time evolution 
with decreasing time steps St = 10 _1 , . . . , 10~ 8 . We perform 
the evolution until the energy is converged up to 12 digits and 
the entanglement entropy up to 6 digits. The run times of the 
algorithm depends very strongly on the model and the param- 
eters used iflOTI . 

We will limit ourselves in this paper to cases in which the 
iTEBD algorithm converges effectively to the matrices of a 
given dimension that are the best approximation of the ground 
state. In a gapped phase, the entanglement entropy of the ex- 
act wavefunction is finite, and the approximation by MPS's 
converges rapidly l30ll . The iTEBD result should therefore 



converge to the physically correct distribution if the correla- 
tion length £ (i.e., the correlation length in the physical sys- 
tem, not the finite-entanglement approximation) satisfies 



a < £ < logx, 



(ID 



where a is the correlation length and x i s the matrix size. 
When these conditions are satisfied, the iTEBD result is not 
evolving rapidly with x, and the numerical spectrum should 
be that of the physical system. The approach of the numeri- 
cal spectrum to the physical one is somewhat subtle, because 
the physical spectrum even at finite correlation length has an 
infinite number of nonzero eigenvalues, while the numerical 
spectrum has at most x- However, the fraction of total proba- 
bility in the tail of small eigenvalues becomes very small with 
increasing x away from criticality, unlike at the critical point. 

Note that this limit is different from the finite-entanglement 
scaling limit log x *C £, where the physical correlation length 
is effectively infinite but there is an apparent correlation length 
in numerics. Hence a different type of result is expected when 
the correlation length satisfies £ ^> logx- It was observed pre- 
viously 12811 that the Calabrese-Lefevre distribution approxi- 
mately describes the spectra obtained in this limit, where the 
correlation length is generated by finite x- In the follow- 
ing section, we compare numerical spectra to the Calabrese- 
Lefevre distribution and in some cases to exact results. 



IV. NUMERICAL RESULTS 

A. Transverse Ising model 

We start by considering the quantum Ising model whose 
Hamiltonian is 



with the Pauli matrices af and g > 0. The model is critical 
at g = 1 with central charge c = 1/2. This model belongs 
to a class for which, the entanglement spectrum can be calcu- 
lated by exploiting the relation between quantum chains and 
two-dimensional classical models i2, HH HH). Using a cor- 
ner transfer matrix approach, it can be show that the reduced 
density matrix has a diagonal form 



p = /Ccxp -^ej-ctc 



(12) 



3>0 



with fermionic operators cj , c- and single particle eigenvalues 

£ , = { (2 ' + 1)e (13) 

where I = 0, 1, 2, The value of e is given by 

e = wl(k')/l(k) (14) 
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where I(k) is the complete elliptic integral of the first kind 
and k' = yl — k 2 . The parameter k, is for the case of the 
transverse Ising model, simply given by 

k={ i/g '- 9> ; as) 

[a ,.9<i- 

We compare the eigenspectrum uj n of the reduced density ma- 
trix p in Eqn. ( fT2l ) with our numerical results. FIG. Ola) shows 
the good agreement of the entanglement spectra obtained from 
Eqn. (TTZt and the numerical results which we obtained using 
the iTEBD algorithm. FIG. |TJb) shows single particle eigen- 
values and compares the analytic results from Eqn. dT2b with 
the numerical results. In particular, we observe a linear disper- 
sion which is predicted by the CTM approach. If we approach 
the critical point, the value of x has to be chosen increasingly 

10° % ■ ■ ■ , 



large to get a good agreement. Our findings agree with previ- 
ous numerical studies for these models in Refs. fToLl25l[32l 

Next we calculate the distribution function n{u>) using the 
iTEBD algorithm and the CTM approach. The results are 
shown in FIG. |2] At the critical point (g=l), the correlation 
length £ diverges and the actual distribution of the eigenvalues 
is flat. The iTEBD algorithm, however, truncates the entangle- 
ment spectrum, thereby generating a finite correlation length. 
Thus the entanglement entropy is finite [28] and we use the 
resulting b = S/2 as a parameter in n(u>). If we detune the 
system from criticality, the correlation length is finite and the 
iTEBD is essentially exact for a x which is large enough. Con- 
sequently, the iTEBD algorithm and the CTM approach give 
the same results. The distribution function n(a) in Eqn. (0 
agrees well with the iTEBD results, independent of whether 
the correlation length is induced by the truncation or an ac- 
tual finite correlation length. If the systems is strongly de- 
tuned from criticality, the discrete structure of the entangle- 
ment spectrum - which is not captured by Eqn. (0 - becomes 
more pronounced. 
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B. XXZ model 



The XXZ spin chain is described by the Hamiltonian 



and is critical in the entire range 7 £ [0, 1] with central charge 
c = 1. Critical exponents change continuously with 7 113311 . 
The results for the distribution function n(u>) near the Heisen- 
berg (XXX) point at 7 = 1 are shown in FIG. [3] For all param- 
eters, the distribution functions Eqn. (0 describes the iTEBD 
results well. A detailed study of the entanglement spectrum 
of XXZ chains can be found in Ref. B4. 
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FIG. 1: Entanglement (Schmidt) spectra of the transverse field Ising 
model for different parameter. Panel (a) compares the numerical 
iTEBD-results (x) with the entanglement spectrum obtained by the 
CTM approach (o) (2I. [3l|l . Panel (b) compares directly the single 
particle eigenvalues. 



FIG. 2: Comparison the prediction by the formula n(A) = 
7o(2\/&log(A max /A)) with the results obtained by the iTEBD al- 
gorithm and the CTM approach for the transverse Ising model near 
the critical point at g = 1. 



5 



10' 




iTEBD (7*1.0, z=100) 
iTEBD (y=0.9, z=100) 
iTEBD (7=0.8, X =10m 
l [2[blog(a> max /a»] 1 ^ 



10" 



1/(0 



FIG. 3: Comparison the prediction by the formula n(oj) — 
Io(2-\/61og(a; ma x/w)) with the results obtained by the iTEBD al- 
gorithm for the XXZ model near the critical point at 7 = 1. 
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C. Spin-1 model 
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FIG. 4: Comparison the prediction by the formula n(u) = 
Io(2^/felog(c<; ma x/u;)) with the results obtained by the iTEBD al- 
gorithm for the spin-1 model near the critical point at 6 — —n/4 



all parameters, the distribution functions Eqn. (O describes 
the iTEBD results well. 



D. Moments of the reduced density matrix 
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FIG. 5: The quantity c a = R a /e- b{a - 1/a) with b = -logcj raax 
calculated numerically for critical points of different models. 

The derivation of the distribution function n(u>) in Eqn. (O 
has been made under the assumption that the moments of the 
reduced density can be expressed as R a = c a e~ b< ' a ~ 1 ^ a ^ 
with c Q being constant. In order to check the validity of 
this assumption, we calculated the moments R a = Tip a ex- 
plicitly using the iTEBD algorithm. The results for c Q = 
iVe- fc(Q ~ 1/Q) and X = 100 are shown in FIG.|5] We find 
c a = 1 in two cases: (i) If a = 1 because Trp = 1 (ii) 
For a — * 00 because in this case, R a = w" lax = e~ ab . In 
between, the c a deviate from unity by at most 10%. This 
explains the good agreements of the iTEBD results and the 
distribution function in Eqn. 



V. SUMMARY 



The two models we have studied so far are integrable mod- 
els. Now we consider the 5=1 Heisenberg chain with bi- 
quadratic term, 

H = ^2 [cos^(S i .S i+1 ) +sm#(S ? ,S 4+1 ) 2 ] , 

i 

which is generally non-integrable. This system is known to 
have two exactly solvable critical points, the SU(2)2 point at 
6 = -tt/4 with c = 3/2 JH and the SU(3)i point at 9 = 
7r/4 with c = 2 Bo] . The entire region 9 e [tt/4, 7t/2) is 
critical with c = 2 07[|. A detailed study of critical SU(iV) 
can be found in Ref. 0811 . The distribution functions n(uS) 
near the critical point at 9 = —ir/4 are shown in FIG.|4] For 



For all considered models, the distribution function n(u>) 
derived by Calabrese and Lefevre agrees at the few percent 
level with the iTEBD results, independent of whether the cor- 
relation length is induced by the truncation or an actual finite 
correlation length. A direct calculation of the moments c a 
suggests that the assumption c a = const, is a good approxi- 
mation. For future work, it would be interesting to understand 
theoretically why this is the case. In general, one would ex- 
pect nonuniversal quantities such as the c a not to have any 
constraint on their variation, in which case it is puzzling why 
the standard models we have studied all have roughly constant 
c Q . Conversely, if there are actually bounds on the c Q that 
explain the observed behavior (that they are not constant but 
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nearly so), then understanding the origin of those bounds is an 
important problem. It is interesting to note in FIG. [5] that the 
model with the strongest deviations from the constant behav- 
ior of moments c a is the spin-1 model, which is the most com- 
plicated analytically and also the most difficult numerically 
(in part because of its higher central charge, which means that 
more degrees of freedom are critical). 

The results of this paper show that the Calabrese-Lefevre 



distribution is indeed a good starting point, at least for a num- 
ber of integrable and non-integrable models, for applications 
of the entanglement spectrum at critical points. They also 
place constraints on the ultimate theory of this entanglement 
spectrum including non-universal effects. The authors ac- 
knowledge support from ARO (F. P.) and NSF DMR-080441 3 
(J. E.M.). 
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